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A  large  number  of  the  Army's  present  and  future  optica1,  svstems 
have  as  their  primary  tasic  the  i^ge  formation  of  a  distant  object. 

In  virtually  all  of  these  systems,  atmospheric  turbulence  plays  a  role 
in  limiting  the  image  quality.  This  is  so  because  the  atmospheric 
medium  is  inhomogeneous  in  temperature  and,  hence,  in  refractive 
index.  Thus,  virtually  all  military  optical  systems  must  operate  in 
media  that  vary  randomly  in  space  and  time.  Even  nonimage-forming 
optical  devices  such  as  laser  rangefinders^  are  affecteo  by  the 
grosser  manifestations  of  curbulence  which  take  the  form  of  erergy 
redistribution^)  wirhin  the  optical  beam  and  raudot.,  pointing  of  the 
beam.1')'  However,  for  device;  such  as  telerccpes  th*t  *re  uaed  to 
form  an  image,  the  random  modu\ation  of  the  optical  phase  as  light 
traverses  the  medium  acts  in  a  more  subtle  way  to  lower  the  resol  v.'ng 
ability  of  ar.  instrument  below  it£  ir,  vacuo  (diffra*  tion-Jiml  ed)  per¬ 
formance  level. 

Because  of  the  deleterious  effect  of  turbulence  on  rhe  perform¬ 
ance  of  Image-forming  instruments,  we  have  reexamined  the  techniques 
of  image  formation  in  terms  of  a  new  approach  to  a  problem  that  has 
afflicted  astronomers  for  over  two  centuries.  We  will  review  the 
basic  mathematical  relationships  central  to  classical  optical  imagery 
along  with  the  constraints  intrinsic  to  the  detection  of  radiation  at 
optical  wavelengths.  It  will  be  shown  that  the  only  parameter  acces¬ 
sible  to  measure  is  the  intensity  (mean  square  cf  the  electric  field) 
and  that  this  parameter  corresponds  to  a  second-order  correlation  of 
the  electric  field. 

Having  developed  the  basis  of  classical  imagery,  we  will  intro¬ 
duce  the  notion  of  higher-order  correlations  of  optical  fields  (fourth, 
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sixth,  etc.)  and  discuss  the  ways  in  which  these  measures  relate  to 
standard  image  formation  with  which  we  are  familiar.  Specifically,  by 
some  manipulation  of  the  far-field  intensity  distribution  of  a  source, 
we  wish  to  infer  the  intensity  distribution  over  the  source  itself. 

We  will  discuss  some  of  the  constraints  of  higher-order  correlation 
techniques,  the  potenrial  they  hold  for  reducing  the  effects  of  trans¬ 
mission  in  the  atmosphere,  and  a  basic  experiment  that  illustrates  the 
mathematical  relationships. 

In  this  paper  we  intend  to  present  a  mathematical  framework 
upon  which  we  can  hang  physical  and  intuitive  arguments.  Although 
extensive  mathematical  treatment  will  be  avoided,  annotation  has  been 
provided  for  those  who  might  wish  to  explore  these  topics  in  detail. 


1.  THE  HUYGENS-FRESNEL  EQUATION 


We  begin  our  discussion  of  optical  imaging  by  examining  the 
coordinate  planes  of  Fi$ .  1.  The  object  to  be  imaged  here  is  found  in 
the  t-n  plane  at  the  leJ t-hand  portion  of  the  figure  and  is  bounded  by 
the  aperture  denoted  I.  In  the  general  case,  the  object  may  be  a 
primary  source  of  light  ->r  a  secondary  surface  from  which  light  is 
scattered.  The  object  is  described  by  a  two-dimensional  source,  since 
an  extension  along  the  axis  of  propagation  can  be  accounted  for,  here, 
by  an  equivalent  field  at  the  plane  of  radiation.  It  is  our  interest 
to  describe  the  nature  of  the  electric  field  as  it  propagates  from  the 
t-n  to  the  x-y  plane,  a  distance,  R0,  away.  For  the  moment,  we  are 
concerned  only  with  the  form  of  the  electric  field  V(x,y)  as  it  relates 
to  the  field  at  the  object,  V(£,n).  The  particular  length  between 
specific  points  in  the  two  planes  is  indicated  by  the  path  R(£,n,x,y). 


Upon  realizing  that  each  point  in  the  aperture  radiates  a 
spherical  wave  to  the  right,  we  can  write  the  Huygens-Fresnel 
principle^  '  for  paraxial  waves  for  which 


V(x,y;t) 


y/‘v(C.r;t) 


expf ikR(s, r;x,y) 1 
iAR(£,n;x,y) 


d£dn. 


(1) 


where  k  =  2*/*,  and  A  equals  the  wavelength  of  the  radiation.  Equa¬ 
tion  (1)  states  tha .  the  electric  field  V(x,y;t)  is  formed  of  a  super¬ 
position  (linear  summation)  of  waves  emanating  from  each  point  within 
the  aperture,  properly  phase-shifted  according  to  the  exponential  term 
and  diluted  by  the  1/R  expression.  Now 

RU,n;x,y)  -  [Ro2  +  (x-o2  +  (y-n)2]**  (2a) 
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,(5) 


when  approximated  by  the  first  two  terms  of  the  binomial^J/  (the  Fres¬ 
nel  approximation)  expansion.  Using  Eq.  (2b)  in  Eq.  (1)  and  making 
the  far-field  (Fraunhofer)  approximation  [R0  >>  k(i’;2  +  n2)  / 2 ]  ,  we 

crn  write  ’  max 

exp  (ik(R  +  (x2  +  y2)/2R  )] 

V(x,y;t)  =  - 


i>.R 


a. 

JJ V(C,  n;t)  exp£- 


t)  exp|-i  —  (x£  +  yn)  |d£dr,. 
Ko 


>]■ 


(3) 


Apart  from  the  coefficients  before  the  integral,  Eq.  (3)  chows  that 
the  electric  fields  in  an  aperture  and  the  far  field  are  related  by  a 
spatial  Fourier  transform  operation.  This  property  is  basic  to  the 
ensuin®  work  involving  intensity  correlations  with  quasimonochromatic 
light.  We  note  that  Eq.  (3)  is  linear;  that  is,  the  two-dimensional 
integral  operator  on  V(£,n;t)  is  linear.  Therefore,  if  either 
V(4,n;t)  or  V(x,y;t)  is  known,  the  other  is  specified  throug.i  a  linear 
(and  hence,  invertible)  transformation. 


2.  THE  VAN  CITTERT-  ZERMIKE  THEOREM 

Although  Eq.  (3)  holds  for  radiation  of  general  frequencies, 
the  relationship  it  expresses  is  somewhat  acade  ic  here  from  the  stand¬ 
point  that  it  describes  electric-field  quantities  that  are  unmeasurable 
at  optical  wavelengths.  No  known  detector  can  follow  oscillations  at 
frequencies  of  1014  hz.  At  optical  frequencies,  the  parameter  acces¬ 
sible  to  measurement  is  the  intensity,  the  mean  square  of  the  electric 
field  averaged  over  many  oscillations.  It  is  this  constraint  that  has 
prompted  many  investigators  to  couch  optical  theory  in  the  form  of 
correlations  of  field  quantities.  Ihe  best  known  of  these  correla¬ 
tions  is  called  the  mutual  coherence  function,  f Qt j , ac2 , t ) ,  where 

/T 

V(x,  ,t  +  t)V*(x2,t)dt,  (4a) 

"  V(Xj ,t  +  f)V*(x2,t) ,  (4b) 

and  t  is  the  time  delay  between  the  instantaneous  product  of  the  elec¬ 
tric  fields  at  the  points  X.J  and  x.„  (and  ^  =  xj.y^).  A  special  cose 
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ot  the  mutual  coherence  function  results  when  the  time  delay  Is  set  to 
zero  or 

1  Oil  *— 2  ' T  ^  t  =  (j  r  ^  CSq  »—2 '  '  (5) 

where  J(x,  ,x,)  is  called  the  mutual  intensity  function. 

a  <L  J 

We  wish  to  calculate  the  mutual  intensity  in  the  far  field  of  a 
spatially  incoherent  source.  Using  Eqs.  (5),  (4b),  and  (3),  we  write 


V(xj ,Vj ;t)V*(x^ ,y, ; t)  =  J(Xj ,y  ;x£ .y,) 

exp[lT  (x!  -  \  •  >i  -  y~->] 


(6a) 


I'R  o)' 


7///  V(fM  ,r,j  ;t)v*(c2,n2;t) 


*  exp  -l 


H  F“  (X1  ’1  "  XA  +  >V'i  ~  v:rl2)|df.1df,2dnIdn2, 


(6b) 


where  the  processes  of  temporal  averaging  and  spatial  integration  have 

been  interchanged  and  the  radiation  assumed  quasimonochromatic 

(-  •  v  _ -  )  so  that  the  wavelength  dependence  can  be  approxi- 

min.  max  min  r  r 

matei  by  the  mean  wavelength.  Since  the  source  is  spatially  incoherent 

the  mutual  intensity  takes  the  forra^) 


V  (  ^  1  >  r ;  >  t )  V*(  C  2  ,  r>2 ;  t)  =  J((.t  ,r'j  ;C2,  r2)  , 

*  -  ^)^(hj  -  r;. 


(7a) 

(7b) 


Physically,  Eq.  (7b)  implies  that  the  time  fluctuations  of  the  elec¬ 
tric  fi  Ids  at  two  non- ident ical  points  in  the  source  plane  are  com¬ 
pletely  uncci related ; *  equivalently,  the  total  power  measuied  at  a 
point  in  the  far  field  is  simply  the  sum  of  the  squared  electric 
fields  from  each  differential  element  of  the  source,  taken  with  the 
proper  phase  delay  and  attenuation. 


If  Eq.  (7b)  is  used  in  Eq.  (6),  the  mutual  intensity  collapses 
to  a  single  aiea  integral,  giving 


The  overbar  of  Eq.  (7a)  by  the  definition  of  Eq.  (4)  indicates  a  time 
average.  Although  spatial  incoherence  is  usually  defined  by  a  time¬ 
averaging  process,  it  can  also  be  defined  in  terms  of  a  spatial 
average  in  the  source  plane. 
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J(Xj  .y^x^y,,)  = 

o; 

yy  1(4, n)  exp|-i  {(Xj-x2)4  +  (y j  -  y7) n]  |d4dn. 


(B) 


Although  usually  presented  in  a  normalized  form,^  Eq.  (8)  is  a  state¬ 
ment  of  the  van  Ci tter t  -  " irnike  theorem  which,  aside  from  the  accom¬ 
panying  coefficient,  shows  that  the  mutual  intensity  in  the  far  field 
of  a  spatially  incoherent  source  is  given  by  the  Fourier  transform  of 
the  intensity  across  chat  source. 


The  mutual  intensity  function  of  Eq.  (£)  is  directly  accessible 
through  measurement  using  o  technique  known  as  Young's  double  pinhole 
experiment.  Figure  1  shows  an  opaque  surface  erected  in  the  x-y  plane. 
Two  pinholes  are  made  in  the  surface.  A  distance  behind  the  surface  a 
fringe  pattern  can  be  observed  oa  the  u-v  plane.  The  fringe  contrast 
relates  to  the  amplitude  of  J(x; ,y j ;x; ,y2> ,  and  the  fringe  shift  cor¬ 
responds  to  the  relative  phase  between  the  radiation  on  the  pinholes. 

If  this  measurement  is  made  for  all  pinhole  spacing*  and  orient  itions 
in  the  x-y  plane,  J(xj ,y i ;x2,y; )  is  completely  specified  and  can  be 
Fourier  transformed  to  derive  I(4,n>,  the  intensity  cn  the  object. 

This  theorem  is  basic  to  most  op.*  (.cal  imaging,  since  lenses  effect  the 
inverse  transformation  of  Eq.  ^4)  to  give  a  scaled  distribution  of  the 
object  i^radiance.  Referring  ogaln  to  Fig.  1,  we  note  that  a  lens  has 
been  placed  its  own  focal  lengt*;  from  both  the  x-y  and  u-v  planes.  We 
assert  without  proof  (see  P^f.  -4,  n.  S3)  that  ir.  this  configuration 
there  is  a  Fourier  transform  rt'.ation  between  these  two  planes  due  to 
the  operation  of  the  lens  on  the  radiation  from  the  x-y  surface  simi¬ 
lar  to  the  relationship  between  the  C-n  and  x-y  planes  located  a  great 
distance  apart.  Thus,  under  the  assumptions  given  in  the  mathematical 
treatment  above,  the  electric  field  V(u,v)  is  t  scaled  distribution  of 
the  field  V(4,r).  The  inversion  of  the  u-v  axes  is  a  mathematical  con¬ 
venience  (see  Ref.  4,  p.  167). 

By  the  scaling  property  from  one  domain  to  another  intrinsic  to 
the  Fourier  transform  operation,  large  detail  in  the  source  is  given 
by  small  separations  in  the  far  field  and  vice  versa.  Thus,  imaging 
with  a  finite  aperture  implies  a  finite  limit  to  the  high-frequency 
detail  resolvable  on  the  source.  This  justifies  the  well-knowr.  de¬ 
scription  of  lenses  as  low-pass  filters.  Also,  since  there  are  many 
more  pairs  of  points  within  an  aperture  corresponding  to  small  separa¬ 
tions  than  to  large,  there  is  a  built-in  redundancy  weighted  in  favor 
of  low-frequency  resolution. 


Equation  (8),  although  a  statement  of  second-order  correlation, 
is  basic  also  to  fourth-order  correlation,  as  we  shall  see  later,, 

3.  THE.  EFFECTS  OF  COHERENCE 


At  thJs  point  in  the  development,  it  is  important  that  we  con¬ 
cern  ourselves  with  the  coherence  properties  of  optical,  sources. (^,9) 

We  again  refer  to  Fig.  1.  Let  us  imagine  that  there  is  a  single  point 
scatterer  at  the  object  (£->*)  plane  illuminated  by  an  optical  source. 
The  field  at  the  x-y  plane  (plane  of  observation)  will  be,  of  course, 
a  spherical  wave.  If  we  begin  to  insert  arbitrarily  other  point  soat- 
tercrs  at  the  object  plane,  the  electric  field  at  the  x-y  plane  will 
be  the  sum  of  all  the  individual  contributors  as  determined  by  geometry 
pec  Eq.  (3).  At  a  particular  instant  of  time,  the  contributions  from 
the  various  scatterers  will  add  constructively  at  some  points  on  the 
x-y  plane,  while  at  other  points  they  may  not.  But,  as  indicated 
earlier,  the  electric  field  is  not  accessible  to  measurement;  the 
intensity  (mean  square  of  the  elect, tic  field)  must  be  averaged  over 
many  rycles  of  field  oscillation.  If  the  optical  bandwidth  of  the 
light  ilium* na ting  the  scatterers  is  sufficiently  narrow  that  the  rela¬ 
tive  phase  across  the  wavefront  remains-,  tons  tone  during  the  period  of 
observation,  then  the  intensity  patter,,  that  would  be  observed  during 
just  a  few  temporal  cycles  of  oscillation  will  remain  flx?d  during  the 
fot.al  period  of  observation-:  Thu  pattern  i&  called  a  speckle  pattern. 
If,  cn  the  other  hand,  the  bandwidth  of  the  source  is  broadened,  then' 
the  speckle  pattern  begins  to  change  with  rime  and  gradually  wasnes 
out  unless  the  period  c£  observ?ticn  ic  shortened  accordingly. 

Fluctuations  of  the  speckle  pattern  in  the  time  domain  relate 
to  the  temporal  coherence  properties  of  t'u-  source,  fluctuations  of 
the  speckle  pattern  in  space  (across  tbo  x-v  plane)  relate  tc  the 
source  sJzs  and  ace  termed  spatial  coherence  effects,  it  is  these 
latter  coherence  properties  that  arc  embodied  in.  the  mutual  intensity 
function  Eq.  (8)  a r.d  relate  directly  to  tt:o  parameter  of  interest, 

I  fjj ,  the  intensify  distribution  over  the  source  plane.  But  the 
ability  to  record  the  speckle  pattern  without  its  being  washed  out 
depends  or-  she  source  bandwidth  *»nd  detection  time  in  the  x-y  plane. 

Speckle  patterns  are  tnsiiy  observed  by  eye  when  laser  radiation 
is  seas  rated  from  a  rough  surface.  However,  with  thermal  sources 
great  must  be  taken  tc-  utilize  c»ea  a  fraction  of  the  speckle  sig¬ 

nal  that  irs  usually  time-variant  with  extremely  short  tLat  constants. 
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4.  A  THEOREM  CONCERNING  THE  FOURTH-ORDER  GAUSSIAN  RANDOM  PROCESS 


We  have  seen  in  Eq.  (7)  the  way  in  which  a  time  average  can  be 
used  to  impose  a  condition  on  the  correlation  of  an  electric-field 
pair.  No  particular  assumption  was  made  about  the  statistics  of  the 
field  variables.  However,  it  is  well  known  that  classical  thermal 
sources  exhibit  statistical  fluctuations  (temporally)  that  are  gaussian 
in  nature.  Hodara'*®'  has  asserted  that  lasers  with  but  a  few  axial 
modes  are,  to  a  good  approximation,  gaussian  as  well.  However, 
Troupe)  has  argued  that  gaussian  statistics  are  achieved  only  in  the 
limit  of  a  large  number  of  axial  modes. 

The  argument  for  gaussian  statistics  can  be  extended  to  the  spa¬ 
tial  domain  as  well.  It  has  been  argued  that^^)  the  received  field 
at  ar.y  point  in  the  far  zone  (as  described  above  in  Section  3)  consists 
of  a  sur  of  random-amplitude,  random-phase,  complex  phasors  .ottributed 
by  the  elementarc  scatterers.  If  the  size  of  the  scattering  area  is 
large  enough  to  include  many  point  scatterers  (or  there  are  enough 
elementary  coherence  areas  composing  'he  source),  the  Central  himit 
Theorem  may  be  used  to  conclude  that  the  electric  field  in  tie  detec¬ 
tion  plane  is  a  gaussian  random  process  in  a  spatial  sense. 


A  well-knowr.  property  of  gaussian  statistics  is  that  all  higher- 
order  moments  are  representable  in  terms  of  the  first  and  second. d^) 

Ir.  particular,  the  fourth-order  correlation  of  electric  fields  (second- 
order  correlation  of  intensities)  car.  be  shown  to  be'1^ 


T  T 


1  T  .  + 


19) 


where  the  defining  relation  of  Eq.  (4a)  has  been  used.  Equation  (9) 
therefore  describes  the  relationship  between  intensity  correlations 
and  field  correlations  for  a  process  that  is  gaussian  in  the  time 
domain.  We  note  that,  in  general,  the  field  correlation  is  a  complex 
quantity,  so  thor  only  the  relationship  between  the  intensity  corre¬ 
lation  and  the  modulus  of  the  field  correlation  is  implied.  It  reveals 
the  underlying  principle  by  which  intensity  correlations  ir.  the  far 
field  may  be  used  to  infer  the  accompanying  field  correlations  and 
hence,  through  Eq.  (8),  to  gain  knowledge  of  the  intensity  distribution 
at  the  source. 


5.  ASSUMPTIONS  BASIC  TO  STATISTICAL  AVERAGING 

We  are  now  in  a  position  to  examine  the  nature  of  the  averaging 
processes  fundamental  to  a  number  of  ootical  processing  schemes. 

Figure  2  represents  ar.  ensemble  of  similar  experiments.  As  in  Section  3, 
we  consider  the  intensity  distribution  in  the  x-plane  to  be  made  over 
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an  observation  time  short  compared  with  any  temporal  fluctuations  in 
intensity.  The  amplitude  distribution  over  the  source  plane  is  identi¬ 
cal  from  one  sample  to  the  next.  The  pail  of  receiver  points  in  the 
detection  plane  also  remains  constant.  However,  each  source  plane  ex¬ 
hibits  a  statistically  similar,  but  independent  phase  mapping  across 
its  extent.  This  insures  that  the  speckle  (intensity)  patterns  observed 
at  the  receiver  (x)  planes  r» main  statistically  identical  as  well.  In 
the  language  of  probability  theory,  each  pair  of  measurements  in  the 
>t-plane  represents  one  particular  sample  in  the  outcome  space  of  the 
experiment . 

We  now  turn  our  attention  to  the  classic  relation  of  imaging, 

Eq.  (8),  It  must  be  recognized  that  there  are  two  fundamental  assump¬ 
tions  basic  to  its  development  which  were  not  explicitly  mentioned. 
First,  the  mutual  coherence  function  of  Eq.  (4)  is  fundamentally  defined 
in  terms  of  an  ensemble  average.  This  is  the  average  that  would  be 
calculated  if  the  total  outcome  space  for  a  given  e/periment  were  known 
for  a  given  space-time  set  of  boundary  conditions.  In  practice,  the 
ensemble  average  is  never  measurable  because  of  a  lack  of  access  to  an 
unlimited  number  of  experimental  configurations  operating  under  identi¬ 
cal  circumstances.  Instead,  the  ergodic  hypothesis  is  invoked:  that 
is,  if  one  particular  experiment  is  performed  under  >.nchan',ing  condi¬ 
tions  (stationarity) ,  then  the  average  that  is  taken  in  the  domain  of 
stationary  conditions  is  assumed  equivalent  to  the  ensemble  (true) 
average.  The  working  definition  of  Eq.  (4)  necessarily  assumes  that 
the  ensemble  average  can  be  replaced  by  a  time  average  because  of  a 
condition  of  stationarity  in  the  time  domain,  the  domain  in  which  the 
averaging  is  made.  Equation  (.8)  is  then  a  time  average  over  many 
instantaneous  products  of  electric  field.  For  this  situation.  Fig.  2 
can  represent  a  time  series  of  field  pairs  (i,j,...n)  which  are  sam¬ 
pled,  multiplied,  and  averaged. 

b,  FOURTH-ORDER  CORRELATIONS  IN  THE  TIME  DOMAIN 

Until  about  two  de.ades  ago,  essentially  all  imaging  was  accom¬ 
plished  by  means  of  second-order  correlations  taken  over  a  time  average 
long  compared  vrith  the  coherence  time  of  the  radiation.  About  that 
rime,  two  astrororers,  Hanbury  Brown  and  Twiss,  were  searching 

for  a  technique  to  infer  the  diameter  of  stars  which  would  not  have  the 
sensitivity  to  atmospheric  turbulence  and  instrument  vibration  charac¬ 
teristic  of  the  stellar  interferometry  of  Michelson  and  Pease.'  J  To 
that  end,  they  were  first  to  suggest  the  use  of  fourth-order  field  cor¬ 
relations.  The  key  to  this  approach  is  embodied  in  Eq.  (9). 

Although  Hanbury  Brown  and  Twiss  did  not  explicitly  use  the 
gauasian  theorem  of  Eq.  (9),  it  is,  in  fact,  intrinsic  to  their 
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mathematical  development.  Their  result  shows,  nevertheless,  that  the 
time-averaged,  two-point  correlation  of  intensities  in  the  far  field 
ot  a  spatially  incoherent  source  is  proportional  to  the  square  of  the 
Fourier  transform  of  t>'°  sluice  intensity  distribution.  We  note  as 
before  that  I'  j  2  is>  *n  general,  a  complex  function  which,  in  fact, 
represents  the  amplitude  spectrum  of  the  source  intensity  distribution. 
Since  in  Eq.  (9)  ' ne  mutual  coherence  function  is  squared,  this  method 
yields  the  power  spectrum  of  the  source  intensity  distribution.  Also, 
since  the  relative  magnitudes  of  the  real  and  imaginary  components  of 
the  amplitude  spectrum  are  unknown,  no  linear  transformation  is  pos¬ 
sible  to  infer  the  intensity  distribution  of  the  source  itself. 


Turning  again  to  Fig.  2,  the  Hanbury  Brown  -  Twiss  experiment  can 
be  understood  by  letting  each  member  of  the  ensemble  represent  a  pair 
of  intensity  measurements  made  in  a  time  less  than  the  coherence  time 
of  the  radiation.  As  many  intensity  products  are  averaged  in  time,  the 
function  approaches  the  ensemble  average. 


Although  this  method  of  imagery  yields  information  only  about 
the  object  intensity  power  spectrum,  it  has  proved  useful  in  the  meas¬ 
urement  of  star  diameters, for  which  only  the  first  zero  crossing 
of  the  Tourier  transform  need  be  known.  Since  the  optical  phase  is 
discarded  immediately  upon  detection,  the  instrumentation  is  not  only 
insensitive  to  mechanical  vibrations  but  to  turbulence- induced  phase 
fluctuations  as  well.  The  primary  limitation  to  the  method  is  the  low 
level  of  radiation  that  is  detected.  As  a  result,  this  approach  is 
limited  to  relatively  intense  stars  to  overcome  the  severe  signal-to- 
noise  problem.  Analyses  of  this  problem  have  been  given  by  Gamo 
and  Twiss. (2i) 


7.  INTENSITY  CORRELATION  IN  THE  SPATIAL  DOMAIN 


With  the  utilization  of  the  laser,  the  experimental  con  traints 
are  quite  different..  The  signal-to-noise  ratio  can  be  increased  typi¬ 
cally  by  six  orders  of  magnitude.  With  a  view  to  exploring  this 

property,  Deitz  and  Carlson^^  have  recently  investigated  the  peten- 
:ial  for  intensity  correlation  in  terrestrial  imagery.  In  many  situ¬ 
ations  of  interest  to  the  \rmy ,  the  availability  of  laser  illumination 
g^eat ’  /  reduces  the  constraints  due  to  noise  limitations.  In  addition, 
since  ground-to-ground  imagery  suffers  from  the  greatest  deterioration 
in  image  quality  because  of  turbulence,  it  offers  the  greatest  poten¬ 
tial  for  improvement. 

(23) 

Although  the  details  of  this  work  can  be  found  elsewhere, 
we  will  again  turn  to  Fig.  2  for  the  basis  of  this  recent  approach.  In 
the  case  of  spatial  intensity  interferometry,  the  signal  in  the  x-plane 
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is  detected  over  a  large  area,  rather  thar  at  just  two  points.  Con¬ 
ceptually,  a  large  piece  of  film  is  placed  in  the  Jc-plane.  The  signal 
is  recorded  spatially  over  a  time  short  compared  with  the  coherence 
time  of  the  radiation  and  is  then  autocorrelated.  This  operation  cor¬ 
responds  to  an  average  taken  in  the  space  domain.  Here,  through  the 
assumption  of  spatial  stationarity  (and  a  sufficiently  large  area  of 
spatial  averaging),  the  spatial  average  is  assumed  equivalent  to  the 
ensemble  average. 

Although  the  method  of  signal  processing  is  very  different  from 
the  Hanbury  Brown  -  Twiss  experiment,  the  final  mathematical  relations 
are  similar.  The  principal  .esult^2^)  shows  that  the  spatial  intensity 
correlation  in  the  far  fie’d  of  a  source  can  be  written 

<I(*|  .T)1(S2.T»  -  j  P-  /  sinc(pT)H(p)dp  (10a) 

0  — 00 
00  u 

«//  expj^-i  “  (Xj-x2)  ‘ijdi 

~oo 

x  fj  C(f)  expj^-i  (Xj  +  x2)  *f.Jdf  | 

—00 

■  I  i  [7  <5,  -S2>]  I  I  S  [£  (i,  +*2>]  |2-  UOb) 

The  angle  brackets  indicate  a  spatial  average,  c  is  the  vacuum  velocity 
cf  light,  r  is  the  average  distance  between  object  and  detection  planes, 
and  w  is  the  angular  frequency  of  radiation.  The  parameter  p  is  the 
difference  frequency  between  modes  in  the  incident  radiation,  the  H 
function  is  defined  in  Ref.  23  and  relates  to  the  temporal  spectrum  of 
the  carrier,  and  k  is  the  wave  number  of  light  (2n/A).  Equation  (10) 
shows  that  the  spatial  intensity  correlation  is  proportional  to  two 
functions.  The  first  is  the  square  of  the  Fourier  transform  of  the 
source  intensity  distribution  as  related  in  Eq.  (9)  by  a  temporal  aver¬ 
aging  process.  The  second  term,  C(f_),  is  a  normalized  phase  corre¬ 
lation  function^12)  describing  the  coherence  interval  over  a  rough 
surface.  In  order  to  describe  a  spatially  incoherent  surface,  C(jf) 
(which  is  defined  by  a  spatial  average)  is  often  allowed  to  assume  the 
role  of  a  delta  function.  In  that  idealized  limit,  its  transform 
becomes  a  constant  and  does  not  band  limit  the  detectable  spatial  fre¬ 
quency  spectrum  of  the  object. '23) 

8.  AN  EXPERIMENT  IN  SPATIAL  INTENSITY  INTERFEROMETRY 

To  illustrate  the  results  of  Eq.  (10),  a  pair  of  crossed  Ronchi 
rulings,  shown  in  Fig.  3,  were  used  as  the  object  in  the  following 
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experiment.  First,  in  order  to  demonstrate  the  Fourier  transform 
relation  of  the  Huygens-Fresnel  equation  [Eq.  (3)],  the  object  of 
Fig.  3  was  transilluminated.  In  the  far  field,  the  intensity  distri¬ 
bution  was  recorded  and  is  shown  in  Fig.  4.  This  signal  corresponds 
to  the  square  of  electric  field  represented  by  Eq.  (3);  equivalently, 
the  signal  of  Fig.  A  is  the  power  spectrum  of  the  electric  field  dis¬ 
tribution  ov^v  the  object  of  Fig.  3. 

To  illustrate  the  results  of  intensity  correlation,  a  seition 
of  ground  glass  was  inserted  (to  gain  spatial  incoherence)  adjacent  to 
the  object  of  Fig.  3  and  again  illuminated  by  a  laser.  The  far-field 
speckle  pattern  that  resulted  is  shown  in  Fig.  5.  The  pattern  appears 
random.  The  mathematical  effect  of  the  ground  glass  has  been  to  multi¬ 
ply  the  object  field,  V( £, n),  inside  the  integral  of  Eq.  (3)  by  a  random 
phase  function. 

Next,  the  speckle  pattern  of  Fig.  5  was  used  to  make  an  identi¬ 
cal  pair  of  optical  transparencies.  These  were  then  used  in  an  optical 
autocorrelator  to  record  the  result  shown  in  Fig.  6.  As  described  by 
the  relation  of  Eq.  (10),  Fig.  6  gives  the  power  spectral  density  of 
the  intensity  distribution  over  the  object.  Since  the  electric  field 
amplitude  and  the  intensity  over  the  object  of  Fig.  3  are  related  by  a 
squaring  process,  the  results  of  Figs.  A  and  6  are  similar  in  the  man¬ 
ner  indicated  by  Eqs.  (3)  and  (10). 

9.  THE  INVERSION  OF  THE  POWER  SPECTRUM 

As  we  have  shown  in  the  previous  sections,  intensity  correla¬ 
tions  lead  only  to  the  power  spectral  density  of  the  source  intensity 
or,  equivalently,  the  modulus  of  the  Fourier  transform.  Depending  upon 
the  application,  the  power  spectral  density  or  its  Fourier  transform, 
the  correlation  function,  may  be  the  parameter  sought  for  system 
use. (24)  jhjg  is  0ften  the  case  in  a  target-guidance  or  recognition 
problem.  However,  without  the  phase  information,  the  Fourier  inversion 
cannot  be  taken  to  derive  the  source  distribution  itself. 

In  the  experiments  of  Hanbury  Brown  and  Twiss,  the  loss  of  phase 
is  not  a  serious  limitation  since  their  objective  is  simply  the  meas¬ 
urement  of  star  diameters.  If  a  circular  disk  is  used  as  a  model  for 
a  star,  the  object  is  known,  a  priori,  to  be  symmetrical.  Thus  the 
spatial  transform  of  the  (real)  intensity  is  pure  real.  For  this  situ¬ 
ation,  the  phase  of  the  transform  is  zero  or  *  for  all  spatial  wave 
numbers,  and  the  square  root  of  the  power  spectrum  can  be  taken  (with 
a  sign  ambiguity)  to  derive  the  spatial  transform  itself. 

£ 

The  author  is  indebted  to  N.  A.  Peppers,  Stanford  Research  Institute, 
for  a  number  of  helpful  suggestions  in  this  experiment. 
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In  the  application  of  intensity  interferometric  techniques  to 
terrestrial  imaging  systems,  though,  the  loss  of  phase  is  a  serious 
limitation  to  the  mechod.  A  number  of  authors  have  addressed  them¬ 
selves  to  the  problem  of  phase  recovery. (25)  /  general  solution  to 

this  problem,  however,  has  never  been  found. 

/  «  / y 

Deitz  and  Carlson  have  recently  proposed  a  scheme''  1  whereby 
the  incoming  electric  field  is  preprocessed  before  intensity  detection 
and  autocorrelation  By  this  method,  the  source  is  effectively  symme¬ 
trized  in  its  intensity  distribution.  If  a  function  is  real  and  ex¬ 
hibits  even  symmetry,  '*7)  its  Fourier  transform  must  be  pure  real. 

Thus,  even  though  the  power  spectrum  is  finally  derived  by  the  measure¬ 
ment  method,  its  accompanying  amplitude  spectrum  is  known  (by  virtue 
of  the  preprocessing)  to  be  pure  real.  Thus,  the  square  root  of  the 
power  spectrum  can  be  taken  to  infer  the  amplitude  spectrum.  The*;*  is 
a  critical  choice  of  signs,  however,  to  be  made  in  the  square  root  pro¬ 
cess.  The  utility  of  the  method  appears  to  rest  with  the  degree  to 
which  system  noise  can  be  precluded  from  interfering  with  those 

decisions. '26) 


10.  CONCLUSIONS 

There  are  a  number  of  special  benefits  from  detecting  images  by 
the  technique  of  intensity  correlation.  (1)  The  method  is  relatively 
insensitive  to  the  effects  of  atmospheric  scintillation. (19)  (2)  Be¬ 

cause  the  signal  is  detected  in  the  spatial-transform  domain,  high- 
frequency  detail  abcut  the  scattering  surface  translates  to  large 
spatial  lags  in  the  far  field.  This  result  could  be  particularly 
important  at  frequencies  where  detector  resolution  is  not  well  devel¬ 
oped.  (3)  A  special  advantage  to  intensity  interferometry  in  the 
spatial  domain  is  the  utilization  of  gaussian  statistics  in  the  spatial 
(not  temporal)  .*ense.  By  this  method,  sources  with  non-gaussian  time 
statistics  (such  as  single-axial-mode  lasers)  can  be  utilized. 

(4)  Still  another  advantage  of  spatial  detection  is  that  images  of 
moving  surfaces  can  be  formed  using  brief  exposures. 

Against  these  benefits  must  be  weighed  the  limitations  of  inten¬ 
sity  interferometry.  The  primary  factor  in  this  respect  rests  upon  the 
greatly  diluted  energy  density  at  the  plane  of  detection.  Because  the 
signals  are  of  such  relatively  low  intensity,  they  tend  to  be  masked  by 
noise.  The  limitation  can  result  from  quantum  noise  in  the  carrier, 
detector  noise,  or  stray  signals  from  unwanted  background.  It  remains 
to  be  seen  whether  these  limitations  can  be  successfully  overcome  so 
that  the  potential  benefits  of  intensity  correlation  can  ultimately  be 
realized. 
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Fig.  1.  Coordinate  axes  for  the  object  (£-n)  and  primary 
detection  vx-y)  planes.  In  addition,  a  lens  is  placed  its 
focal  length  behind  the  x-y  plane.  Finally,  a  focal  length 
behind  the  lens  is  the  u-v  plane. 


Fig.  2.  An  ensemble  of  similar  experiments.  Source 
amplitude  and  receiver  points  are  identical.  Random 
phase  variations  over  each  source  plane  are  statistically 
independent. 
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Fig.  3.  Optical  Transparency 
Used  as  Object 


Fig.  4.  Spatial  Power  Spectral 
Density  of  the  Electric  Field 
Distribution  of  Fig.  3 


Fig.  5.  Far-Field  Speckle  Pattern 
of  Transparency  in  Fig.  3  When 
Used  with  Random  Phase  Screen 


Fig.  6.  Spatial  Power  Spectral 
Density  of  the  Intensity  Dis¬ 
tribution  of  Fig.  3 
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